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RESEARCH MEMORANDUM 


AERODYNAMIC CHARACTERISTICS AT TRANSONIC AND SUPERSONIC 
SPEEDS OF A ROCKET-PROPELLED AIRPLANE CONFIGURATION 
HAVING A 52.59 DELTA WING AND A LOW, 

SWEPT HORIZONTAL TATL 


By Alan B. Kehlet 
SUMMARY 


A flight investigation over a Mach number range from 0.79 to 1.83 
has been conducted in order to determine the aerodynamic characteristics 
at low lift of a rocket model of an airplane configuration having а 
52.59 delta wing of aspect ratio 3.08 with NACA 65A005 airfoil sections 
in the streamwise direction and a low, swept horizontal tail. The 
lift-curve slopes and static longitudinal stability were nonlinear with 
lift coefficient over most of the Mach number range and increased with 
increasing lift coefficient over the lift range covered. The minimum 
drag coefficient increased from 0.016 to 0.035 through the transonic 
speed range. 


The model damping characteristics were irregular and altered for 
the two tail settings over the Mach number range covered with the higher- 
lift tail setting having the greater damping. 


The measured periods of the lateral oscillations were of the same „+ 
order of magnitude as the longitudinal periods, but apparently were not 
affected by lift coefficient. The model exhibited stable static direc- 
tional characteristics throughout the Mach number range tested. 


INTRODUCTION 


As part of a general research program investigating longitudinal 
stability of wings having various plan forms and thickmess ratios, a 
rocket-propelled model of an airplane configuration having a 52.5° delta 
wing of aspect ratio 3.08 has been flown. The basic fuselage-empennage 
configuration (ref. 1) had swept horizontal and vertical tails with the 
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all-movable horizontal tail mounted in a low position. During the 
flight, the horizontal tail was deflected in a square-wave program 
between 0.100 and -3.2°. 


The model was flown at the Langley Pilotless Aircraft Research 
Station at Wallops Island, Va. 


SYMBOLS 
CN normal-force coefficient, эп 8 | 
Сс chord-force coefficient, ES 
CL lift coefficient, CN cos a - Cc sin a 
Ср drag coefficient, Co cos а + Cy sina 
Ст lift coefficient at minimum drag coefficient . 
Cm | pitching-moment coefficient about center of gravity 
an normal acceleration determlned from accelerometer, ft/sec? 
84 longitudinal acceleration determined from accelerometer, 
ft/sec — 
g acceleration of gravity, ft/sec? 
q ` dynamic pressure, O.'TOpM? _ КШ = 
p free-stream static pressure, lb/sq ft 
M Mach number = 
R Reynolds number, based on wing mean aerodynamic chord 
8 wing area (including area enclosed by fuselage), sq ft 
c wing mean aerodynamic chord (M.A.C.), ft 
b wing span, ft | — 
W weight, Tb I Е 
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A cross-sectional area, sq in. 

r equivalent radius of cross-sectional area, in. VA/x 

X longitudinal distance from station O, in. 

1 length of model, in. 

P period of oscillation, sec 

11/2 time to damp to one-half amplitude, вес 

a angle of attack, deg 

Š control panel deflection (measured free stream normal to 


wing chord plane), deg 


Ө amgle of pitch, radians 
B angle of sideslip, deg 
d 
Cm, = 8 
а a( | 
(S 
dCm 
Cmo = : = 
(Fi 
Саз» effective rate of change of yawing-moment coefficient with 
sideslip angle per degree (derived as in ref. T) 


аст Е аст 
2 
аСт, а (ст : — 
Subscripts: 
w wing 
£ fuselage 
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. Lb de i 
x — == (=. 
57.3 dt 
_ dê : 
q = at 


The symbols a, б, а, 5, and B% used as subscripts indicate / 
the derivative of the quantity with respect to the subscript; for 


= 
example, vhs = ^ " — — = 


MODEL AND INSTRUMENTATION 


Model 


A three-view drawing of the model is shown in figure 1. The non- 
dimensional equivalent body and area distribution, presented for 


possible drag correlation purposes, are shown in figure 2. Photographs 
of the model are shown in figure 5. : 


The empennage has a vertical fin of wood and aluminum with the 
quarter-chord line swept 60? апа NACA 654005 airfoil sections in the 
streamwise direction; a horizontal tail of duralumin with the quarter- . 
chord line swept 45° with 20° negative dihedral and NACA 654006 airfoil 
sections in the streamwise direction perpendicular to the chord plane. 
The fuselage is described in reference 2. 


The steel delta wing of aspect ratio 3.08 had a leading-edge sweep 
of 52.59 and NACA 65A003 airfoil sections in the streamwise direction. 


Each panel of the horizontal tail was deflected in an approximate 
square-wave program by & separate servo-control fed by a common pressure 
System and regulated by an electric motor-driven selector valve. For 
the present investigation, the stop positions were 0.1° and -3. 29 


measured parallel to the model center line amd normal to the wing chord 
plane. 


The model weighed 118.75 pounds and had a moment of inertia in 
pitch and yaw of 8.21 and 8.31 slug-ft*, respectively. The center of 
gravity was located at 0.26 of the wing mean aerodynamic chord. 


vel yki, 1 dı1 
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Instrumentation 


The model was equipped. with an NACA telemetering system which 
transmitted continuous measurements of normal acceleration at the 
center of gravity, normal acceleration at a reference nose station, 
angle of attack, longitudinal acceleration, transverse acceleration, 
control position, total pressure, and reference static pressure. 


Flight-path information was obtained from tracking radar and 
atmospheric conditions at altitude from a radiosonde released immediately 
after the flight. Motion-picture cameras were used to photograph the 
launching and first portion of the flight. 


TESTS AND DATA REDUCTION 
Preflight Tests 


Prior to instrumentation, the model was placed in a profile machine 
and measurements were made in the streamwise direction of the airfoil 
sections, perpendicular to the chord planes of the vertical and hori- 
zontal tails and of the wing. These measurements were used to check 
construction tolerances. 


Prior to flight testing and with the instruments installed, the 
model was suspended by shock chords and vibrated by an electromagnetic 
shaker. The following model natural frequencies and modes of vibration 
were determined: 


Second 
bending, 
cps 


Torsion, 
cps 


Horizontal tail . 
Vertical tail. . 
Wing . +. 4 « « « 


Wing Modes 


First bending Second bending Torsion 
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Flight Tests 


The model was launched at an angle of approximately 609 from the 
horizontal by means of a mobile launcher as shown in figure 3. Two 
6-inch-diameter solid-fuel ABL Deacon rocket motors boosted the model 
to maximum velocity. All measurements used were taken during the 


decelerating portion of the flight. - | — | — 


Data Reduction 


The response of the model to deflections of an all-movable hori- 
zontal tail in an approximate square-wave program was analyzed by the 
method of reference 5. The indicated angles of attack were corrected 
to angles of attack at the model center of gravity by the method of 
reference 4. The two-accelerometer method for obtaining instantaneous 
total pitching-moment coefficients was used as described in reference 2. 


ACCURACY = 


The absolute accuracy of the measured quantities is impossible to 
establish because the instrument calibrations cam not be checked during 
or after the flight. Most of the probable instrumentation errors occur 
as errors in absolute magnitude. Incremental values or slopes should, 
in general, be more accurate than the absolute values. An indication 
of the systematic instrument errors possible is given by the followlng 
table, based on an accuracy of £1 percent of the full instrument range: 


£0 .00490 +0.00122 
+.02251 + .00563 
+.04180 +.01045 


The CW Doppler radar unit is believed. to be accurate to better 
than 1 precent for nonmaneuvering models. The Mach number at peak 
velocity should, therefore, be accurate to l. percent or better. Since 
the Mach number subsequent to peak velocity was determined from 
telemetric data, it probably becomes less accurate as the Mach number 
decreases. : — 


Further errors in the aerodynamic coefficients may arise from 
possible dynamic-pressure inaccuracies which are approximately twice 
as great as errors in Mach number. _ 


Sabes EM 


CU 
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Am indication of random errors encountered may be noted from the 
scatter of data points shown in the figures. Errors in angle of attack 
and control panel deflection are independent of dynamic pressure and 
are not likely to vary with Mach number. The horizontal-tail deflec- 
tions are estimated to be accurate within +O. 109 and the increments in 
angle of attack within t0.209. 


RESULTS AND DISCUSSION 


Dynamic pressure amd Reynolds number obtained during the flight 
are shown in figures 4 and 5, respectively. A typical time history at 
the higher Mach numbers of some of the quantities measured is shown in 
figure 6. Hereafter; S in order to distinguish the two tail settings, 
the deflection of 0.1° shall be referred to as the low-lift tail setting 
and the deflection of -3.20 as the higher-lift tall setting. 


Longitudinal Trim 


The variations of the trim lift coefficient and trim angle of attack 
at the two control settings as functions of Mach number are shown in fig- 
ure T. At transonic speeds amd with increasing Mach number, the model 
exhibited a trim change of approximately i nose-up in the low lift 
range; whereas at the higher trimmed lift condition, а nose-down trim 
change of approximately 29 occurred. 


Lift 


The variation of the lift coefficient with angle of attack over the 
Mach nunber range is shown in figure 8. The lift-curve slopes repre- 
sented by the faired lines in figure 8 are presented as functions of 
Mach number in figure 9 for both tail settings. Included in figure 9 
for comparative purposes are the lift-curve-slope dat& for & wind-tunnel 
model with a similar wing-fuselage combination having no horizontal tail 
(ref. 5). In view of the differences in models, the agreement between 
the rocket-model data and the wind-tunnel data is considered good. It 
is indicated in figure 9 that the lift coefficient does not vary line- 
arly with the angle of attack, particularly at the lower Mach numbers; 
however, reference to the trim lift coefficients (fig. 7), which are 
indicative of the lift ranges covered by the two tail settings, shows 
that the linearity of the lift-curve slopes present at the higher Mach 
numbers may be due to the low lift range of the higher-lift tail setting 
and not entirely to a Mach number effect. The increase in lifi-curve 
slope with increasing lift coefficient, where nonlinearity occurs, is 
believed to be due in part to a decrease in the downwash field over the 
low horizontal tail with increasing angle of attack of the model and to 
а nonlinearity in the lift-curve slope of the long nose section of the 


fuselage. 
ee 
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Drag К 


The variation of drag coefficient with lift coefficient corresponding 


to the lift ranges of figure 8 is shown in figure 10. The maximum lift- 
drag ratios that could be measured from figure 10 and lift coefficient 

at which (L/D)max occurs are shown in figure 11. Because of the —__ 
reduced amplitude of the oscillations, the model never reached its maxi- 
mum lift-drag ratios at Mach numbers greater than 1.12; values of 
(L/D)max and Cr, for (L/D),,, obtained at higher Mach numbers were _ 
extrapolated from the data for the drag due to lift, minimum drag coef- 
ficient, and lift coefficient at minimum drag. At high subsonic speeds, 
the maximum value of the maximum lift-drag ratios obtained was about 9, 
decreasing to about 6 at low supersonic speeds and to about hz at the 


test's limit; the values of lift coefficient corresponding to these 
values are about 0.25, 0.40, and 0.30, respectively. 


The minimum drag coefficient and the lift coefficient at. minimum 
drag coefficient obtained from figure 10 are presented as a function 
of Mach number in figure 12. The flagged symbols on tbe minimum-drag- 
coefficient curve indicate extrapolated points from the higher-lift 
tail setting. It should be noted that, as indicated by the basic data 
of figure 10, the values of the Êh Mi arig coefficient correspond 
closely to the values of the zero lift-drag coefficient, the drag polars 
at the low-lift tall setting are flat around zero lift coefficient, and 
the values of the lift coefficient at minimum drag can be determined, 
with any degree of accuracy, only below a Mach number of 1.5. At all 
Mach numbers for the higher-lift tail setting, no lift coefficient ai 
minimum drag can be obtained. 


By using the values of the lift coefficient at minimum drag coeffi- 
cient (fig. 12) for the low-lift tail setting and by interpolating 
values for the higher-lift tail setting, the effect of lift on drag ав 
а function of Mach number was determined and is presented in figure 15. 
As indicated by the data of figure 15, the model exhibited poor leading- 
edge suction over the Mach number range where the effect of lift on drag 


could be determined. As stated before, prior to instrumentation, measure- ` 


ments of the airfoil sections were made. The wing airfoil sections, as 
determined from these measurements, revealed a leading-edge radius some- 
what smaller then the true NACA 65A005 airfoil section. The resulting 
rather sharp leading edge is believed to be par: of the. cause of the 
lack of leading-edge suction obtained. 


Longitudinal Static Stability 


The measured periods of the longitudinal oscillations for both 
tail settings as a function of Mach number and the lóngitudinal-static- 
Stability parameter Сто determined from these periods are presented in 


cour; 
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figures lh and 15, respectively. As with the lift-curve slopes (fig. 9), 
the nonlinearity present in the data decreased with increasing Mach num- 
ber. Where nonlinearity of the data occurred, the higher-lift tail 
setting had the greater value of static stability.  Agsin, however, the 
linearity at the higher Mach numbers may be due to the decreasing trim 
lift coefficient rather than to the increasing Mach number. 


The variation of the static-stability parameter is reflected in 
the aerodynamic-center location for the two tail settings (fig. 16). 
The data at both tail settings exhibited & rearward movement of the 
aerodynamic center up to & Mach number of about 1.5 with the higher- 
lift tail setting having the greater stability. Above M- 1.5, & 
forward movement occurs. 


The variation of the total pitching-moment coefficient, obtained 
from the two-accelerometer method (ref. 2), with lift coefficient is 
shown in figure 17. Although some scatter is present, particularly at 
subsonic speeds, the data agree, in general, with the slopes indicated by 
the period method and the trim data of figure T which are also shown 
in figure 17. 


A measure of the horizontal-tail effectiveness in producing moment 
and ability to produce lift, as obtained from an average of Cm, for 


the two tail settings (fig. 15), the trim angle of attack (fig. T), апа 
an estimated tail length, is presented in figure 10. Both parameters 
exhibit the same general shape as the lift-curve slopes (fig. 9), that 
is, increasing with Mach number to M = 1.0 and then decreasing with 
increasing Mach number. At the test limit (M = 1.85), & decrease of 
approximately 50 percent of that at М = 1.0 is noted in both param- 
eters. Included in figure 18 for comparative purposes are Су and 
Cms, Obtained from Cry, from the exposed duralumin wing of reference 2 
(extrapolated Cr, at the higher Mach numbers) and corrected for 
dihedral effects. 


Although the values obtained from reference 2 are approximately 
20 percent higher throughout the Mach number range covered, it should 
be pointed out that the duralumin wing CL used was not corrected 
for the small gap that exists between the fuselage and horizontal-tail 
panel. The effects of a small gap in a nonviscous flow (ref. 6) can 
be large, reducing the lift-curve slope by as much as 40 percent. By 
using an average gap of 0.076 inch and Cr, from reference 2, an esti- 


mate of the effect of the gap at M = 1.00 was made. The gap reduced 
the parameters by about 30 percent (about 15 percent lower than the 
model values). However, since the flow is viscous and, therefore, the 
effective gap is smaller than the measured gap, the tail parameters 
obtained from Cre of reference 2 and corrected for an effective gap 


should, in general, be close to the values obtained from the model. 


i RET LL 
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Damping in Pitch ЕГ | B 
The variations of the amplitude ratio (the ratio of the amplitude j 
from the trim line of successive oscillations to the initial oscilla- 
tion) with time, for both tail settings, are shown in figure 19. The 
time to damp to one-half amplitude, represented by the time required 
for the faired line to cross 0.50 amplitude ratio (fig. 19), as a 
function of Mach number, amd the pitch damping-moment factors are shown 
in figure 20. Ag with the lift-curve slope and period data, nonlinearity 
with the trim lift coefficient was also present in the time to damp to 
one-half amplitude over all but the higher Mach numbers. Where non- 
linearity occurred, the model exhibited irregular damping character- 
istics with the higher-lift tail setting having the greater damping. 
Since, for the configuration used in this investigation, damping is _ 
primarily due to the tail and would be expected to decrease as downwash 
over the tail decreases, and, since, іп the discussion of lift, downwash 
was believed to decrease over the tail with increasing angle of attack, 
the model damping at the higher-lift tail setting would be expected to . 
be less than tbe damping at the low-lift tail setting. The reasons for 
the altered damping characteristics are not known at this time; however, 
since the model exhibited both longitudinal and lateral oscillations, 
the possibility exists of coupling between these oscillations. The 
effect of coupling is believed to tend to reduce the mddel damping. 


Directional Static Stability ° 


Although the primary purpose of the model flown in this investiga- 
tion was longituđinal stability, lateral oscillations were present and, 
therefore, were investigated. These oscillations are believed to be 
caused by some asymmetry in the horizontal-tail stops shat could be 
present amd cause lateral disturbances. 


The variation of the measured periods of the lateral oscillations 
as a function of Mach number and the static-directional-stability coeffi- 
cient (derived as in ref. T) as determined from these periods are pre- 
sented in figures 21 and 22, respectively. Throughout the Mach number 
range the measured periods of the lateral oscillation are of the same 
order of magnitude as the measured periods of the longitudinal oscilla- 
tions and appear to be unaffected by lift coefficient over the lift 
range covered. Although some coupling between the longitudinal and 
lateral oscillations is probably present, the effect of this coupling 
is believed to be small on all the aerodynamic parameters with the 
exception of the damping in pitch. 


The static-directional-stability coefficient exhibited a decrease 
of about 30 percent from the high subsonic to the highest supersonic 
Mach numbers. Flexibility becomes increasingly important with increasing 


«онаш . — š 
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Mach number, and the decrease in the static-directlonal-stabllity coef- 
ficient is believed to be due in part to flexibility of the vertical 
tail. Reference to the table of natural frequencies indicates that the 
vertical tail is rather flexible. At the Mach numbers covered, however, 
the model exhibited stable static directional characteristics. 


An estimation of the approximate maximum amplitudes of the angle 
of sideslip over the Mach number range covered has been made. At the 
high subsonic Mach numbers, maximum amplitudes of +20 and 150 occurred 
for the higher- and low-lift tail settings, respectively; these ampli- 
tudes decreased to about half at the transonic and supersonic Mach 
numbers. 


CONCLUSIONS 


A flight investigation of the aerodynamic characteristics at 
transonic and supersonic speeds of a rocket-propelled airplane con- 
figuration having a 52.59 delta wing and a low, swept horizontal tail 
indicated the following conclusions: 


1. At transonic speeds and with increasing Mach number, the model 
exhibited a trim change of approximately 1° nose-up in the low lift 
range and а 2° nose-down trim change while trimmed at the higher lift 
condition. 


2. The lift-curve slopes were nonlinear in the transonic and low 
supersonic Mach number range and increased with lift coefficient over 
the lift range covered. 


3. The minimum drag coefficient increased from 0.016 to 0.035 
through the transonic speed range. The minimum drag coefficient was 
close to the drag coefficient at zero lift. 


lh. Over the Mach number range where it could, be determined, the 
model exhibited poor leading-edge suction. 


5. The aerodynamic center moved rearward with increasing Mach num- 
ber up to a Mach number of about 1.5, with greater stability at the 
higher lifts. 


6. The model damping characteristics were irregular and altered 
for the two tail settings over the Mach number range covered with the 
higher-lift tail setting having the greater damping. 
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7. The model exhibited stable static directional characteristics 
over the Mach number and lift ranges covered. 


Langley Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., January 7, 1954. 
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Figure l.- General arrangement of model. 


All dimensions 


Wing 
Aspect ratio 3.08 
Area (total) 2.976 sq ft 
Area (exp) | .940sq ft 
Dihedral OO de 
B; M.A.C. 131 ft 
Airfoil section МАСА 
14.00 


Horizontal tail 
Aspect ratio 400 
Area (total) 0.905 sq ft 
Dihedral -20.0 deg 
M.A.C. 0.504ft 
Airfoil section NACA 
65A 006 


Vertical tail 
Area (total) | 1.37 sqft 
Airfoil section NACA 
65A 003 


are in inches. 
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Model 


Equivalent body 


$2 90.25 tn. 


Horizontal tail 
Vertical tall 
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Area distribution 


Nondimensional equivalent body and area distribution of model. 
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(8) Three-quarter front view, 19578 61 


Figure 3, Photographs of models. 
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(b) Model on launcher. 


Figure 5.- Concluded. 
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Mach number 


Figure L.- Variation of dynamic pressure with Mach number. 
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Figure 5.- Variation of Reynolds number with Mach number. 
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Mach number 


(a) Trim Lift coefficient. 
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(b) Trim angle of attack. 
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Figure 7.- Longitudinal trim characteristics as в function of Mach number. 
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(b) Supersonic Mach numbers. 


Figure 8.- Variation of lift coefficient with angle of attack at transonic 
and supersonic Mach numbers. 
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Figure 9.- Variation of lift-curve slope with Mach number. 
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Lift coefficient 


Figure 10.- Variation of drag with Lift. 
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(a) Maximum lift-drag ratios. 


Mach number 


(b) Lift coefficients at which maximum lift-drag ratios occur. 


Figure 11.- Variation of maximum lift-drag ratios and lift coefficients 
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at Which maximum lift-drag ratios occur as a function of Mach number. 
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(а) Minimum drag coefficient. 
(b) Lift coefficient at minimm. drag coefficient. 


Figure 12.~ Variation of minimum drag coefficient and lift coefficient 
at minimum drag coefficient as a function of Mach number. 
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Figure 15.- Effect of lift on drag as a function of Mach number. 
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Figure 14.- Variation of measured period with Mach number. 
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Figure 15.- Variation of static-stability parameter with Mach number. 
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Figure 16.- Variation of aerodynamic center with Mach number. 
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Lift coeffioient 


(b) 8 = -5,20, 


Figure 17.- Variation of total pitching-moment coefficient with lift 


coefficient. 
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(a) Ability to produce lift. 


(p) Effectiveness in 


Mach number 


Figure 19.- Variation of ability of horizontal tel 
oducing moment as 4 function of Mach number. 
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Amplitude ratio 
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Time from initial peak, sec 
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(b) 850.19. 


Figure 19.- Variation of amplitude ratio with time. 
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Mach number 


(a) Time to damp to one-half amplitude. 


Mech numoer 


(b) Pitch damping-moment factor. 


Figure 20.- Variation of damping characteristics of short period oscilla- 


tions with Mach number. 
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Figure 21.- Variation of measured period of lateral oscillation with 
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Mach number. 


Figure 22.- Variation of static-directional-stability coefficient with 


